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Abstract 

The Global Colour Model of QCD is used in conjunction with a pure-gluon 
lattice correlator (by Marenzoni et al.) to extract from meson data a 
momentum-dependent quark-gluon coupling down to s « 0.3GeV 2 .This is 
compared with a lattice calculation (by Skullerud) of the quark-gluon coupling. 
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A key feature of Quantum Chromodynamics (QCD) is that the quark-gluon coupling varies 
strongly with gluon momentum q over the range < q < 2GeV relevant to low energy hadronic 
physics. Here we extract from meson data this quark-gluon coupling g{q 2 ) down to q = 0.5GeV, 
using the Global Colour Model (GCM) of QCD in conjunction with a pure-gluon lattice cor- 
relator by Marenzoni et al. 0]. The extracted quark-gluon coupling (see Fig. 2.) is compared 
with a recent lattice calculation by Skullerud |2| . 

The GCM modelling of QCD is based on the idea that because the hadronic correlators 
are given by explicit functional integrals it should be possible, after an appropriate change of 
variables of integration, to identify a dominant configuration. It turns out that the existence of 
this dominant configuration is nothing more than the constituent quark effect. The GCM has 
been recently reviewed by Tandy [|| , and applied to a number of problems [Q, ||, ||, [?], [| |Io|, |ll| . 

In the functional integral approach correlators are defined by 

_ JVqVqVADCVC^.gix) exp(-S Q CD[A,q,q,C,C\) 

JVqVqVAVCVCexp(-S Q cD[A,q,q,C,C}) 

The various complete correlators Q lead to experimental observables. They are related by an 
infinite set of coupled Dyson-Schwinger Equations (DSE), and by the Slavnov- Taylor gauge- 
symmetry-related identities (Non-scripted G"s will denote constituent correlators, as defined 
later). However the GCM does not derive from these equations/identities, its nature follows 
instead from an analytical continuum estimation procedure for the functional integrations. 
Direct numerical estimation procedures are used in lattice modellings of the functional integrals. 
The correlators in (|j]) may be extracted from the generating functional of QCD 

Zqcd[v,V,J] = / VqVqVAVCVCexp(~S QC D[A,q,q,C,C} +m + m + JA). (2) 

The functional transformations which lead to the GCM are discussed in Tandy fl; briefly 
and not showing source terms for convenience, the gluon and ghost integrations are formally 
performed 

J VqDqVAVCVCexp(-S Q cD[A,q,q,C, C]) 
= J VqVqcxp{- J q(-j.d + M)q+ 

+ f / 3% (y)G»* ( X ~y) + §J H3%Q a X + ) (3) 

where j^(x) = q(x)^-j l _ l q(x), go is the bare coupling constant, and Q^ v {x) is the gluon correlator 
with no quark loops but including ghosts 

_ _ fVAVCVCA°(x)A:(y)eM-SQCD[A,C,C]) 
V) ~ JVAVCVCeM-S Q CD[A,C,C]) ' W 

A variety of techniques for computing Q fJl ^(x) exist: the gluon-ghost DSE [Q, and the gluon 
only DSE |l3| and lattice simulations [l4| . The terms of higher order than the term quartic in 
the quark fields are difficult to explicitly retain in any analysis. However we can model, in part, 
the effect of these higher order terms by replacing the coupling constant go by a momentum 
dependent quark-gluon coupling g(s), and neglecting terms like Q^ p and higher order. This g(s) 
is a restricted form of vertex function. This modification goQ filJ (p) — > D^ il ,{p) — g{p 2 ) 2 Q^(jp) 
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and truncation in (JJj) then defines the GCM. However we make one further modification: we 
shall use lattice results with ghosts neglected for Q llv {x) (l). Then g(s) models as well the effect 
of the ghosts in both the gluon correlator and the quark-gluon vertex. See jl2j for an analysis 
of these ghost effects. We call D /liy (p) the effective gluon correlator. 

The GCM is equivalent to using a quark-gluon field theory with the action 



S GCM [A, q, q]= J (q(-j.8 + M+ iA a J^^)q + ^D^(id)A a v 



(5) 



Here D~^(p) is the matrix inverse of D^ v (p), which in turn is the Fourier transform of D fiu (x). 
This action is invariant under q — > Uq,q — > qU\ and A^X a — > UA®\ a U^ (where U is a 
global 3x3 unitary colour matrix) - the global colour symmetry of the GCM. The gluon self- 
interactions that arise as a consequence of the local colour symmetry in (H) and the ghost 
and vertex effects lead to D~^(p) being non-quadratic. Hence, in effect, the GCM models the 
QCD local gluonic action J [A]F* [A], having local colour symmetry, in Sqcd of (|l|), by 
a highly nonlocal action, having global colour symmetry, in the last term of (||). The success 
of this modelling has been amply demonstrated |j. The form for g(p 2 ) is here determined by 
comparing the meson data determined D^ v (p) to the pure-gluon lattice-determined correlator 

Hadronisation of the functional integrations in ([!]) involves a sequence of changes of vari- 
ables involving, in part, the transformation to bilocal boson fields, and then to the usual local 
hadron fields (sources not shown): 

Z « J VqVqVAcxp(-S GC M[A, q, q]+rjq + qiq) (GCM) 
= J DB...exp(-S[B,..}) (bilocal fields) (6) 

= J T>n'Dp'Duj...exp(—Shad['K,p,oj....]) (local fields) . (7) 

The bilocal fields in naturally arise and correspond to the fact that, for instance, mesons 
are extended states. This bosonisation/hadronisation arises by functional integral calculus 
changes of variables that are induced by generalized Fierz transformations that emerge from 
the colour, spin and flavour structure of QCD p5| . The final functional integrations in (^) over 
the hadrons give the hadronic observables, and amounts to dressing each hadron by, mainly, 
lighter mesons. The basic insight is that the quark-gluon dynamics, in (0), is fluctuation 
dominated, whereas the hadronic functional integrations in (Q) are not. 

The second key idea in the GCM is that in proceeding from (||) to ([?]) one expands S[B, ..] 
about the configuration Bcq that minimises it; giving the GCM Constituent Quark (CQ) 
equations. 

S S 

w&yj 1bcq = ° ' (8) 

Thus for all hadrons one assumes a universal dominant configuration. This amounts to assum- 
ing that all hadrons share a common dynamical feature. Of the set B(x,y)cQ only A(x — y) 
and B(x — y) are non-zero translation-invariant bilocal fields characterising the dominant con- 
figuration. Then writing out the translation invariant CQ equations we find that the dominant 
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configuration is indeed simply the constituent quark effect as they may be written in the form 

i, 

G- 1 {p)=ik + ™ + \ J ^D^ip-q^Giq)^, (9) 

and we see that this is the gluon dressing of a constituent quark; and is exact in the GCM. Here 

G(q) = (iA(q)q. 1 + B{q) + m)" 1 = -iq."f(J v (q) + <J s {q). (10) 

In the chiral limit there are more Bqq fields that are non-zero, and a resultant degeneracy of 
the dominant configuration is responsible for the masslessness of the pion |3| . 

The constituent quark G correlator should not be confused with the complete quark corre- 
lator Q from (Q). This Q would be needed to analyse the existence or otherwise of free quarks. 
The G on the other hand relates exclusively to the internal structure of hadrons, and to the 
fact that this appears to be dominated by the constituent quark effect. The evaluation of Q is 
a very difficult task, even in the GCM. G is however reasonably easy to study using <^). 

The hadronic effective action in (^) arises when S[B, ..] is expanded about the dominant 
CQ configuration: the 1st derivative is zero by (||), and the 2nd derivatives, or curvatures, give 
the constituent or core meson correlators G(q,p; P) 

G ' 1 (^ ;P ^ F - r -( ^,^(,,, ) |g -)' (U) 
after exploiting the translation invariance and Fourier transforming. Higher order derivatives 
lead to couplings between the meson cores. The G(q,p;P) are given by ladder-type correlator 
equations, see j|. The non-ladder effects can be inserted by the final functional integrals in (Q), 
giving the complete GCM meson correlators Q(q 7 p;P). In the present analysis the u> and ai 
mesons are described by these constituent meson correlators; that is, we ignore meson dressings 
of these mesons. The mass M of these states is determined by finding the pole position of 
G(q,p;P) in the meson momentum P 2 = —M 2 , this leads to the homogeneous vertex equation 

r(p; p ) = -\J - p^ g (i + f p Wi - f ( 12 ) 

To solve (|) for various D fJiU (p) and then to proceed to use A(s) and B(s) in meson correlator 
equations for fitting observables to meson data is particularly difficult. A robust numerical 
technique is to use a separable expansion for Z? M „ (p — q) . In Landau gauge 

JV(p) = OW - ^ W). and ^W-tV-^fM? 2 )- (is) 

We expand D(p — q) in (^) into 0(4) hyperspherical harmonics 

D{p -q) = D Q (p 2 , q 2 ) + q.pD^p 2 , q 2 ) + ... (14) 

where 

2 /" r 

Do(p 2 ,q 2 ) = - df3sin 2 (3D(p 2 + q 2 - 2pqcos/3), ... (15) 
* Jo 

We then introduce multi-rank separable expansions for each term 

D„( P 2 ,q 2 )= Y, IMp 2 )IM 9 2 ),.... (16) 
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Introduction of the separable expansion clearly breaks translational invariance and must be 
regarded purely as a numerical procedure, much like a lattice breaks translation invariance. 
Translation invariance is restored as the rank of the separability is increased. Here we use a 
rank n = 3 form for Dq, and rank 1 form for D\. The constituent quark equations then have 
solutions of the form 

B{s) = Y d B i {s); B i (s)=b i T i (s),... (17) 
where the 6^, .. are easily determined, in the chiral limit, to be 

6 s 2 = y7r 2 / sdsBi(s)a 8 (s). (18) 

where 
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scr v (s) 2 + a s (s) 
and er s and a v are seen to have the form of sums 



CT s(s) = X! CT s( s )^ &v(s) = ^ a v(s)i, (20) 

i— l,n i—l,k 

However rather than specifying 1^ in ( [l6| ) we proceed by parametrising forms for the a s and 
<r v ; then the T; follow from (0) and (|l9|): 

cr s (s), = Ciexp(-dis), i = 1, 2; cr s (s) 3 = c 3 



2s - <f 3 (l - exp(-2s/d 3 )) 



2s 2 



<r„(s) = —J . (21) 



2s-/3 2 (l-exp(-2s//3 2 )) 

As these forms are entire functions we avoid spurious singularities developing in G. The asymp- 
totic form of (7s (s) ~ 1/s 2 for s — * oo is described by the cr s (s)3 term. With these parametrised 
forms we can numerically relate the mass of the a\ and ui mesons, from ([l2]) and / ff (for Nf = 2) 
to the chiral- limit parameter set {ci,C2,c 3 , , d\, d,2, d%, /3} in a robust and stable manner. The 
parameter values are shown in Table 1. The chiral limit expression for f v is, sec Q, 

u=6 Jwi i ^ ' 2{(Tsa ' s + saWv) ~ s{asa " ~ {a ' s)2) " s2{av< ~ K)2) B{q)2 ' (22) 

The translation invariant form for the effective gluon correlator is easily reconstructed by 
using D(p 2 ) = D (p 2 ,0) from (|J|) 

D{p2) = ? W ^(o? pV„(p 2 )2 + ^(p 2 ) 2 ' (23) 



With the parameter set in Table 1, @ gives &i = 0.0210 Gel/ 2 , 6 2 = 0.0251 GeV 2 and 
63 = 0.0351 GeV 2 and the resulting D(p 2 ) is shown in Fig.l; it has estimated uncertainties 
of 5%. Shown in Fig.l for the pure gluon correlator is T>{p 2 ) from the lattice calculations, 
corresponding to the value j3 = 6.0, of Marenzoni et al (lj where the errors arise from a 5% 
uncertainty in the lattice spacing; a = 0.50 ± 0.025Gel/ _1 . In Fig. 2 we show the form of 
g(s), where g 2 (s) — D{s)/T>{s), see dig)), that then follows from our analysis. Here the error 
bars now indicate combined uncertainties. This extracted quark-gluon coupling extends down 
to 0.3Gey 2 , and shows infrared (IR) enhancement. Below this limit the separable expansion 
becomes unreliable unless more terms and more fitting data are used. We have not corrected 
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for either lattice spacing dependence or for quark loops; corrections for these would require 
further development. It is possible to identify where the IR effect arises. If we artificially lessen 
this effect at small s then we find that the main consequence is an increase in the value of f^. 
Indirectly, then, we can show that the IR signature is the (inverse) pion size in comparison with 
the a\ and li masses. The pion size enters through because in (|2^) in the chiral limit the pion 
form factor T„(q;0) = B(q), see Tandy [[|. However the GCM extraction of this effect does 
not explain what aspect of QCD drives it. In it is argued that the QCD origin of this IR 
enhancement is due to the ghost correlator presence in the quark-gluon vertices. We also report 
various condensate values that arise from the present work: < qq > = (211.4MeV) 3 \iGeV, 
< gqF^a^q >= (491.5MeV) 5 \ 1GeV , and < ^F^F^ >= 0.026GeV 4 \ 1GeV , ignoring quark- 
loop contributions. 

Our most significant result follows from comparing, in Fig. 2, the GCM-meson-data/lattice- 
gluon determined quark-gluon coupling with that determined recently by Skullerud |2| using a 
lattice calculation with (3 — 6.0(a = 0.5GeV _1 ) and a lattice size of 16 3 x 48. For comparison 
we also show the perturbative coupling derived from the two-loop beta function 

g 2 (s) = (& ln(^) + ^lnln(^)) > (24) 

with b = 11/16tt 2 ,6i = 102/(16tt 2 ) 2 for A = 0.420GeV. Fig.2 indicates a general agreement of 
all three methods down to s w 0.7GeV 2 . The most significant difference being the decreasing 
lattice g(s) in the deep IR; however this could be due to the finite lattice size which induces 
an IR cutoff, or to the absence of the ghost effects Q. The Skullerud data is similar to 
the running coupling extracted from the 3-gluon vertex |Q. In Fig. 3 we show a s = g 2 /(47r) 
against q(GeV). These results indicate that QCD may now be sufficiently well modelled by 
the GCM in the low energy regime that detailed hadronic calculations may be performed, 
particularly for the nucleon properties; the GCM having the advantage of easily dealing with 
the near chiral limit needed for the nucleon, in contrast to lattice studies. Fig.2 shows that 
the lattice results for the gluon correlator and the quark-gluon coupling may be combined to 
form D(s)i at — g 2 (s)i a t'D(s)iat] a lattice derived effective gluon correlator for except for 
the deep IR where we should be guided by the meson data fitting. We thus have a meeting of 
the continuum and lattice approaches. The deep IR behaviour remains undetermined, but the 
region of uncertainty mainly affects questions of absolute confinement and will have little effect 
upon low energy hadronic phenomena. 

We thank N. Stella for assistance with the lattice results in [[lj. Research supported by an 
ARC Grant from Flinders University. This work is part of the activities of the Special Research 
Centre for the Subatomic Structure of Matter, University of Adelaide. 
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Table 1: cr s (s) and cr v (s) Parameters 



Cl 


0.1732G6V- 1 


di 


1.389GeV- 2 


C2 


1.527G6V- 1 




4.982GeV- 2 


C3 


0.3435GeV 3 


d 3 


1.971GeV 2 




0.4807GeV 







7 



Figure Captions 



Figure 1 The effective gluon correlator D(s) (solid line) extracted by fitting the GCM to meson 
data. Also shown are the lattice results for the pure gluon correlator V(s) from Marenzoni et 
al. (1995). The error bars indicate uncertainties arising from the value of the lattice spacing 
a = 0.50±0.025Gey- 1 . 

Figure 2 The GCM quark-gluon coupling g(s) (boxes). Here g 2 (s) was obtained by dividing 
the GCM effective gluon correlator D(s) (solid line in Fig.l) by the lattice gluon correlator 
T>{s). The error bars arise from the lattice spacing uncertainty and from systematic errors in 
the fitting of the GCM to the meson data. Also shown is g(s) from the lattice calculation of 
Skullerud (circles) (1997). The curve shows the two-loop form for A = 0.420GeV. 

Figure 3 Here we replot, for convenience, the GCM quark-gluon coupling in the form a s = 
g 2 /{Air) against q(GeV). 
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